Entanglement entropy and multifractality at localization transitions 
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The von Neumann entanglement entropy is a useful measure to characterize a quantum phase 
transition. We investigate the nonanalyticity of this entropy at disorder-dominated quantum phase 
transitions in noninteracting electronic systems. At these critical points, the von Neumann entropy 
is determined by the single particle wave function intensity, which exhibits complex scale invariant 
fluctuations. We find that the concept of multifractality is naturally suited for studying von Neu- 
mann entropy of the critical wave functions. Our numerical simulations of the three dimensional 
Anderson localization transition and the integer quantum Hall plateau transition show that the 
entanglement at these transitions is well described using multifractal analysis. 
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I. INTRODUCTION 

Entanglement is a unique feature of a quantum system 
and entanglement entropy, denned through the von Neu- 
mann entropy (vNE) measure, is one of the most widely 
used quantitative measures of entanglement ,ii2i2i£ Con- 
sider a composite system that can be partitioned into two 
subsystems A and B. The vNE of either of the subsys- 
tems is sa = — T^aPa hi pa = sb = —TtbPb hi pb- Here, 
the reduced density matrix pa is obtained by tracing over 
the degrees of freedom in B: pa = b\iP ab) {ifr ab\ and 
similarly for pg. In general, for a pure state \iPab) of a 
composite system, the reduced density matrix is a mix- 
ture, and the corresponding entropy is a good measure 
of entanglement. 

The scaling behavior of entanglement entropy is a par- 
ticularly useful characterization near a quantum phase 
transition 2 -. The entanglement entropy can show non- 
analyticity at the phase transition even when the ground 
state energy (the quantum analog of the classical free en- 
ergy) is analytic. While these ideas have been studied in 
a number of translation-invariant models,—^ there have 
been far fewer investigations of random quantum critical 
points (for notable exceptions, see Ref. [g). 

In particular, noninteracting electrons moving in a dis- 
ordered potential can undergo continuous quantum phase 
transitions between an extended metallic and a localized 
insulating state as the Fermi energy is varied across a crit- 
ical energy Ec- Well known examples are the Anderson 
transition in three dimensions and the integer quantum 
Hall (IQH) plateau transition in two dimensions where 
the ground state energy docs not exhibit any nonanalytic- 
ity. In contrast, vNE will be shown to exhibit nonanalyt- 
icity at these transitions and a scaling behavior. At the 
outset, it should be emphasized that because of the single 
particle and disorder-dominated nature of these quantum 
phase transitions, entanglement as characterized by vNE 
and its critical scaling behavior are fundamentally differ- 
ent from those calculated for interacting systems. This 
statement will be made more precise later. 

In a noninteracting electronic system close to a dis- 
ordered critical point, the wave function intensity at en- 



ergy E, \tpE(r)\ 2 , fluctuates strongly at each spatial point 
r and, consequently, has a broad (non-Gaussian) distri- 
bution even in the thermodynamic limit.— This non-self- 
averaging nature of the wave function intensity is char- 
acterized through the scaling of its moments. In partic- 
ular, moments of normalized wave function intensity, P q 
(called the generalized inverse participation ratios), obey 
the finite-size scaling ansatz, 



P q (E) = J2 l^(r)| a * ~ L-^ T q [{E E C )L^} . (1) 

r 

Here, L is the system size, v is the exponent characteriz- 
ing the divergence of correlation length, £e ~ \E— Ec\~ v ■ 
T q is called the multifractal spectrum, and the over- 
bar denotes averaging over different disorder realizations. 
T q [x) is a scaling function with Tq(x — ► 0) = 1 close to 
the critical point E = Ec- When E is tuned away from 
Ec, the system either tends towards an ideal metallic 
state with P q (E) ~ L~D( q -i) ^ D being the num b er f 

spatial dimensions) or becomes localized with P q {E) in- 
dependent of L. 

Below, we first show that the disorder-averaged vNE 
can be expressed as a derivative of P q and thus, its scaling 
behavior follows from multifractal analysis. After that, 
we apply our formalism to understand the numerical re- 
sults on vNE at the three dimensional Anderson localiza- 
tion and IQH plateau transitions. vNE in the Anderson 
localization problem was studied previously*^ but the 
connection with mulitfractality and the unique features 
of vNE at these quantum phase transitions have not been 
clearly elucidated. 



II. ENTANGLEMENT ENTROPY IN 
DISORDERED NONINTERACTING 
ELECTRONIC SYSTEMS 

Even though the disorder induced localization problem 
can be studied in a single particle quantum mechanics 
language, there is no obvious way to define entanglement 
entropy in this picture. However (see Ref. 0), entangle- 
ment can be defined using the site occupation number 
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basis in the second-quantized Fock space. Let us divide 
the lattice of linear size L into two regions, A and B. 
A single particle eigenstate of a lattice Hamiltonian at 
energy E is represented in the site occupation number 
basis as 



(2) 



Here iPe{t) is the normalized single particle wave func- 
tion at site r and \n) r denotes a state having n particles 
at site r. We decompose the above sum over lattice sites 
r into the mutually orthogonal terms, 



|V-e) = |1)a®|0) b + |0) a ®|1)b 

where 

|l>x = $>*(r)|l>r®|0> P ,, 



(3) 



A = Q<) \o) r (4) 



with analogous expressions for the \1)b and |0)b states. 
Notice that these states have the normalization 

(0\0) A = (0\0) B = 1, (1\1) A = p A , (1\1) B =p B , (5) 

where 

p A = J2\Mr)\ 2 , (6) 

and similarly for p b with p A + ps = 1 • 

To obtain the reduced density matrix p A , we trace 
out the Hilbcrt space over B in the density matrix 
p = \?Pe)(4>e\- This gives, 



p A = \l) A (l\+p B \0) A (0\. 
The corresponding vNE is given by 

s A = -p A In pa - Pb lnps 



(7) 



(8) 



In the above equation, we see that manifestly s A = sb- 
More importantly, s A is bounded between and In 2 for 
any eigenstate. This is in sharp contrast to the entan- 
glement entropy in interacting quantum systems where 
it can be arbitrarily large near the critical point. The 
reason for this is also clear: Even though we used a 
second-quantized language, we are dealing with a single 
particle state rather than a many body correlated state. 
Consequently, the entanglement entropy does not grow 
arbitrarily large as a function of the size of A. 

We also observe that at criticality, if the whole system 
size becomes very large in comparison with the subsys- 
tem A, we can restrict the subsystem to be a single lat- 
tice site and study the scaling dependence with respect to 
the overall system size L. Then, using the ansatz of scale 
invariance, we can always find the scaling of the entan- 
glement as a function of the subsystem size I since near 
criticality, only the dimensionless ratio L/l can enter any 
physical quantity. To extract scaling, we find the bipar- 
tite entanglement of a single site r with the rest of the 



system and sum this over all lattice sites in the system. 
Using Eq. ((8]), we write this as 



S(E) = -J2{H'E(r)\ 2 ln\i, E (r)\ 2 

reL d 

+ [l-|Mr)| 2 ]ln[l-NMr)| 2 ]}. (9) 



To leading order, the second term inside the square 
bracket in Eq. ([9]) can be dropped since \ipE(r)\ 2 <C 1 at 
all points r when the states are close to the critical energy. 
We can readily relate the disorder average (denoted by 
overbar) of this entropy to the multifractal scaling in Eq. 
([T]) and get the L scaling as 



S(E) 



dPq 

dq 



dq 



oq 



(10) 



We do not know the general form of the scaling function 
T q , but we can get the approximate L dependence of the 
entropy in various limiting cases. For the exactly critical 
case when T q = 1 for all values of q, we get 



S(E) ~ ailnL, 



(11) 



where the constant a\ = dr q / 'dq\ q —i is unique for each 
universality class. From the discussion following Eq. fl}, 
the leading scaling behavior of S(E) in the ideal metal- 
lic and localized states is given by DhiL and ailn^E, 
respectively. From the limiting cases, we see that, in 
general, S(E) has the approximate form 



S(E) ~ fC[(E - E C )L 1 ' U ] lni, 



(12) 



where the coefficient function JC(x) decreases from D in 
the metallic state to ai at criticality and then drops to 
zero for the localized state. We will see that this scaling 
form is verified in our numerical simulations. 



III. ENTANGLEMENT IN THE THREE 
DIMENSIONAL ANDERSON MODEL 

The scaling form for the entanglement entropy aver- 
aged over all eigenstatcs of the single particle Hamilto- 
nian is also of interest since this scaling can change as 
a function of disorder strength. To be specific, let us 
consider the 3D Anderson model on a cubic lattice. The 
Hamiltonian is 



(13) 



where cj(ci) is the fermionic creation (annihilation) op- 
erator at site i of the lattice, and the second sum is 
over all nearest neighbors. We set t = 1, and the Vi 
are random variables uniformly distributed in the range 
[-W/2, W/2]. It is knowr>i£ that as W is decreased from 
a very high value, extended states appear at the band 
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FIG. 1: (Color online) Scaling curve in the 3D Anderson 
model. With the choice of v — 1.57 and C = 12.96, all data 
collapse to a universal functions f±(x). The two branches 
correspond to w < and w > 0. 
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FIG. 2: (Color online) S(E, W, L) as a function of E and W 
computed in a system with L = 10. The square shows the 
mobility edge reported in Ref . 0. Because of the finiteness of 
the system, the transition from the localized to the delocalized 
region is smooth. 



center below the critical disorder strength W c — 16.3, 
and a recent workii reported the localization length ex- 
ponent v = 1.57 ±0.03. 

The analysis leading to Eq. (fT2|) also holds when we 
study wave functions at a single energy, say E = and 
increase the disorder strength in the Anderson model 
across the critical value W c . In this case, the states at 
E = evolve continuously from fully metallic to critical 
and then finally localized, resulting in the approximate 
form for the entanglement entropy as 

S(E = 0,w,L) ~C(wL^ v )]nL, (14) 

where w = (W — W c )/W c is the normalized relative disor- 
der strength and C(x) is a scaling function. In particular, 
as mentioned before, C(x) — > D as w — > — 1, C(x) — > as 
w — > oo, and C(x) = a± when w = 0. 

Next, we look at the energy-averaged entropy. We av- 
erage Eq. pop over the entire band of energy eigenvalues 
and construct the vNE, 

E 

where L 3 is also the total number of states in the band. 
Then using Eqs. (fT2|) and (fT4|) , one can show that close 
to w = 0, 

S(w,L) ^C + L- 1/v f ± (wL 1 ^)\nL, (16) 

where C is an L independent constant, and f±(x) are 
two universal functions corresponding to the two regimes 
w > and w < 0. 

We numerically diagonalizc the Hamiltonian [Eq. [T3] 
in a finite L x L x L system with periodic boundary con- 
ditions. The maximum system size is L = 13, and the 



results are averaged over 20 disorder realizations. The 
scaling form of S(w,L) is given by Eq. (fTfJ|) . Figure [1] 
shows the results of the data collapse with a choice of 
v = 1.57, and the nonuniversal constant C = 12.96 is 
determined by a powerful algorithm described in Ref. fl2l . 
The successful data collapse reflects the nonanalyticity 
of the von Neumann entropy and accuracy of the multi- 
fractal analysis. 

We also use the transfer matrix method^ to study the 
energy dependence of S(E, w, L) by considering a quasi- 
onc-dimcnsional (quasi-lD) system with a size of (mL) x 
L x L, m > 1. Wc use L up to 18, and m = 2000 > 1 
is found to be sufficient. To compute vNE, we divide the 
quasi-lD system into m cubes labeled by I = 1, 2, . . . , m, 
each containing L 3 sites. We normalize the wave function 
within each cube and compute the vNE, S 1 (E,W, L), 
in the I th cube, and finally S(E, W, L) is obtained by 
averaging over all cubes. 

A typical ~S(E, W, L) with L = 10 is shown in Fig. H 
The value of S(E,W,L) is normalized by ln(L 3 ) such 
that S — ► 1 in a fully extended state. The energy E is 
normalized by (W/2 + 6), which is the energy range of 
nonzero density of statesJ^ The mobility edge computed 
in Ref. [13 is also plotted in Fig. O The validity of the 
scaling form in Eq. (fl"4)) is seen in Fig. [3] In particular, 
the function C(x) shows the expected behavior. 



IV. ENTANGLEMENT IN THE INTEGER 
QUANTUM HALL SYSTEM 

Consider now the second example, the integer quantum 
Hall system in a magnetic field B. The Hamiltonian can 
be expressed^ in terms of the matrix elements of the 
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FIG. 3: (Color online) The quantity C in Eq. g4j). The range 
of the system sizes is too small to observe the weak L depen- 
dence. Inset: S(E = 0, W, L) as a function of In L for three 
different W . 



states \n, fc), where n is the Landau level index and k 
is the wave vector in the y direction. Focussing on the 
lowest Landau level n = 0, with the impurity distribution 
V(r)V(r') = V 2 <5(r-r'), the matrix element (0, k\V\0, k') 
can be generated as in Ref. [l(| 

Now, consider a two dimensional square with a lin- 
ear dimension L = \/2~kMIb, where lg = (h/eB) 1 / 2 is 
the magnetic length and M is an integer, with periodic 
boundary conditions imposed in both directions. We dis- 
cretize the system with a mesh of size y/Trl B /V2M. The 
Hamiltonian matrix is diagonalized and a set of eigen- 
states {\ip a } = Ylk a k.a\0, k)}^ =1 is obtained with cor- 
responding eigenvalues {Ea}^^. The energies are mea- 
sured relative to the center of the lowest Landau bancU^ 
in units of T = 2Vq /^/2ttIb- Finally, for each eigenstate 
the wave function in real space can be constructed as 

i>a{x, y) = (x, y\lp a ) = ^2 a h,a1p0,k(x, y), (17) 



where ?po t k{x, y) is the lowest Landau level wave function 
with a momentum quantum number fc. 

The dimension of the Hamiltonian matrix increases as 
Nk ~ M 2 , making it difficult to diagonalizc fully. In- 
stead, we compute only those states \ip a ) whose energies 
lie in a small window A around a preset value E, i.e. 
E a e [E - A/2,E + A/2]. We ensure that A is suf- 
ficiently small (0.01) while at the same time, there are 
enough states in the interval A (at least 100 eigenstates). 

We now uniformly break up the L x L square into 
nonovcrlapping squares Ai of size I X I, where I = 
Ib\Ztt/2, independent of the system size L. For each 
of the states, we compute the coarse grained quantity 
I( X y)eA \^a(x,y)\ 2 dxdy. The computation of the vNE 
for a given eigenstate follows the same procedure de- 



scribed for the Anderson localization. Finally, by av 
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FIG. 4: (Color online) Scaling of the von Neumann entropy 
S(E) for the IQHE. M instead of L is used in the data collapse 
with the accepted value of v = 2.33. 



craging over states in the interval A, the vNE S(E, L) 
is obtained at the preset energy E. The scaling form 
o_f ~S(E, L) is given by Eq. H2J) with E c = and is 
S(E,L) = JC(\E\L 1 /")lnL. A good agreement with the 
numerical simulations is seen in Fig. 2] 



V. CONCLUSIONS 

We have clearly established the formalism for com- 
puting the entanglement entropy near quantum critical 
points in noninteracting disordered electronic systems. 
We have also identified its relation with the well-studied 
notion of multifractality and illustrated our concepts 
through numerical simulations of two important models, 
the 3D Anderson transition and the IQH plateau tran- 
sition. This work represents a starting point to study 
entanglement in electronic systems with both disorder 
and interactions. 
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